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Using the lattice NRQCD action for heavy quark, we calculate the heavy quark expansion parameters fi^ and 
for heavy-light mesons and heavy-light-light baryons. The results are compared with the mass differences 
among heavy hadrons to test the validity of HQET relations on the lattice. 
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1. Introduction 

In the calculation of inclusive decay rates of the 
heavy hadron, the heavy quark expansion (HQE) 
technique is widely used. At the order 1/toq of 
HQE two nonperturbative parameters 



2M 



2Mhc 



He 



He 



Q{iDfQ 



Qa-BQ 



He 



He 



(2) 



appear in the calculation. Hq represents a heavy- 
light meson or heavy-light-light baryon (for h 
hadrons, H^ = B, B* , Af,, Eh, E^). For instance, 
the lifetime ratio of b hadrons is given as H 
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with cg — 1-2. While Hq may be evaluated from 
experimental values of hyperfine splitting, the de- 
termination of /i^ requires some theoretical in- 
puts. It should be noted that the parameters are 
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defined in the static limit: mg — > oo. For heavy- 
light meson, /^^ q has been calculated using the 
lattice version of the Heavy Quark Effective The- 
ory [§. 

In this work we calculate /i^ and /i^ on the 
lattice using the NRQCD action for heavy quark. 
Although the individual matrix element suffers 
from large perturbative uncertainty due to power 
divergence in the matching calculation, their dif- 
ferences like fil{Hl^^)~^j.l{Hl^^) are free from the 
uncertainty of the operator. We calculate both 
^,liHi'^) - f^UHi'^) and ^^UH^'^) - ^,1{H^^\ 
and compare them with the corresponding pre- 
dictions for mass splittings. 

2. HQET mass formula 

The parameters /i^ and //^ can be indirectly 
obtained from hadron masses, using 




where A is the residual energy difference between 
Mhq and mg surviving in the infinite heavy 
quark limit, /ij and /ig. appear in the correc- 
tion terms of 0(l/mg). Therefore, by consider- 
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Figure 1. Hyperfine splitting of ground state 
mesons. Circles is obtained from the energy dif- 
ferences, while crosses are from the matrix ele- 
ments. 

ing proper mass differences, certain combinations 
of /ij and /ig can be extracted. 

For example, a difference of can be obtained 
from the mass splitting in a spin multiplet, be- 
cause A and /i^ have the same value. Also, the 
spin averaged mass — {Mb + 3Af^)/4 does 
not depend on /x^., because is proportional to 
the spin of the light degrees of the freedom and 
the sum of /Iq in the spin multiplet vanishes. 

3. Lattice calculations 

We carry out quenched QCD simulations 



at /3=6.0 



20^ 



48 lattice. 



The 



NRQCD action including all 0(1/toq) terms and 
the non-perturbatively improved clover action 
(csu,=l-769) is adapted for heavy quark and light 
quark, respectively. Five heavy quark masses 
amQ=\.i, 2.1, 3.0, 5.0, and 10.0 are used to study 
the 1 /mg dependence of hadron masses and ma- 
trix elements, while three hopping parameters 
J^=0. 13331, 0.13384, and 0.13432 are simulated 
to extrapolate to the chiral limit i<L'c=0. 135284(8). 
The inverse lattice spacing a~^=1.85(5) GeV is 
determined with the p meson mass mp=770 MeV. 

We measure the three-point functions 
{OHQ{t)O^Mto)Olj^{0)), where Ohq is an 
interpolating field to create or annihilate the 
hadron Hq, and Ott.g is the operator to be mea- 
sured, Q{iD)^Q or Qa ■ BQ. We divide them by 
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Figure 2. Hyperfine splitting of heavy-light-light 
baryons. 

{OHQ{ti)0^u^{Q)) to obtain the desired matrix 
elements /i^ and /ig. 

4. Hyperfine splittings 

From (|^) the hyperfine splitting Mb* — Mb is 
given by — A/iQ/2mQ, or equivalcntly 

Ml. - Ml = EE -{iil{B*) - filiB)), (5) 

at the leading order. In Figure we plot our 
results for —A^q together with the measurement 
of M|. - M|. We observe that the relation (|) 
is satisfied very well, while both are significantly 
lower than the experimental values for B and D 
mesons. 

In deriving (^) we used a relation 

Af,l=t4{B*)-f4iB)=0, (6) 

which holds in the static limit. However, for the 
NRQCD action including the spin-magnetic inter- 
action term at 0{l/mQ), the operator mixes 
with Og at order as/mg. This is the reason why 
our result for — A/i^ deviates from that of the 
mass difference in the lighter heavy quark mass 
region. In other words, the relation (^) may be 
considered as a renormalization condition for the 
operator 0^. 

Similar analysis can be made for the hyperfine 
splitting of heavy-light-light baryon, i.e. E* — 
S splitting. Figure || shows the mass difference 
and the matrix element —A^q. Both arc in good 
agreement. 
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Figure 3. l/A^s dependence of Ma^ — Mq. 



5. Ma, - Mb 

The heavy-light meson-baryon mass difference 
il/A(, — Mjj is given as 



A{Ab) - A(B) 
1 



2mc 



f^liA,)+f^liB)]. (7) 



The intercept at 1/Mg=0 yields A(Ab) - A{B) 
while the slope is described by — A/iJ = 
-{^^liA,)-^^UB)). 

In Figure y we plot Ma, — as a func- 
tion of 1/Mg. For the intercept we obtain 
A(A6) - A(B)=393(31) MeV. in agreement with 
a previous work by Ali Khan et al, A(A;,) — 
A(B)=415(156) MeV. Our resuh is shghtly larger 
than the experimental values for b and c hadrons. 
However, to draw a definite conclusion we have to 
consider several systematic errors, especially the 
finite volume effect, because our lattice may not 
be large enough for baryons. 

The slope obtained from the fit of the mass dif- 
ference is consistent with zero: —0.21(21) GeV^. 
Our results of direct measurement of — A/i^ is 
plotted in Figure ^, which is consistent with the 
result from mass difference, but have much better 
accuracy. Our result is also compatible with the 
phenomenological estimate —0.01(3) GcV^ |^ ob- 
tained from a combination {M\^ — Mq) — (A/ac ~ 
Md)- 
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Figure 4. — A/i^ measured from the matrix ele- 
ments is compared with the slope of mass differ- 
ence A/a, - Mg, -0.21(21) GeV^. 



6. Conclusions 

We confirm that the lattice measurements of 
the matrix elements /z^ and are consistent 
with the HQET mass relations. The well-known 
problem of quenched lattice calculation that the 
hyperfine splitting is much smaller than the ex- 
periments is also reproduced. 

An important extension of our work is to mea- 
sure the matrix elements of four-quark operators, 
which are relevant to the 1 /rriQ corrections to the 
lifetime ratios ||l|. 
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